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Abstract. - A novel pattern forming instability in a mixture of a granular material and water 
in a horizontal rotating drum is experimentally investigated. The particles accumulate in radial 
symmetric rings separated by pure water. The transition between the homogeneous state and 
the structured state is hysteretic. The transition point is extrapolated from the growth rates 
of the chain of rings. The trajectory of a single particle is discussed to estimate an upper 
boundary for the transition point. 



Introduction. - Granular media show an interesting and peculiar dynamic behavior p^, 
and are also of major importance in technological processes. One of the standard problems 
of the powder technology is the mixing of materials jz] . It is often accompanied by pattern 
formation like stratification .'! or segregation 4, 5, 6,, 5 . The axial and radial segregation in 
rotating-drum mixers with dry multicomponent mixtures, like small and large particles or 
particles of different densities, has been investigated. The different angles of repose of the two 
species are considered as the segregation mechanism for the axial segregation. The separation 
of large and small particles has recently been observed in a mixture of particles dispersed 
in a lower density fluid [8]. In the presence of a fluid the time for the pattern formation is 
significantly decreased. 

Granular media immersed in a fluid show a type of segregation even for monodisperse 
particles The volume fraction of granular particles in these experiments is much 

less than in the previous cases. Furthermore, the driving frequency of the drum is sufficiently 
high for the centrifugal forces to become dominant. The observed pattern is a chain of cir- 
cumferential rings of high particle concentration. Tirumkudulu et al. reported on this type 
of axial segregation in partially filled horizontal cylinders In both cases they used 

suspensions of particles, which were neutrally buoyant. If the device was completely filled 
with the suspension, no instability was observed at all. The experiment was modeled theoret- 
ically by Govindarajan et al. ^3] including shear-induced diffusion of particles, concentration- 
dependent viscosity and the existence of a free surface. Boote and Thomas |10|. and Thomas 
et al. |12j reported similar experiments with particle densities less and more dense than the 
surrounding fluid in a partially fluid-filled cylinder. For particles of a different density than 
the fluid, the free surface is not necessary for the instability. Lipson observed this phenomenon 
during crystal growth in a supersaturated solution |14| . 
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Fig. 1 - (a) Experimental setup, (b) Photograph of the tube with the pattern. A top view of the 
tube is shown. The sand is accumulated in the white areas, (c) Corresponding transmitted light 
intensity, sand accumulations appearing dark. The 5 pixel wide frame indicates the area considered 
for further analysis. 



In this publication we present experiments on axial segregation of monodisperse particles 
which have a higher density than the surrounding fluid. The experiments have been carried 
out in a drum, which was completely filled with a sand-water mixture. We estimate an upper 
boundary for the critical frequency for the onset of the pattern formation by the numerical 
analysis of the trajectory of a single particle. 

Experiments. The experimental setup is shown schematically in Fig. ^a). The ex- 
periment consists of a horizontally aligned plexiglass tube, which is 400 mm long. Its inner 
diameter measures 23.73 ± 0.05 mm. Both ends are closed with plexiglass plugs. The tube is 
filled with 2.5% (by volume) glass beads and 97.5 % water. This is sufficient to reach about 
85% of the close packing for a monolayer on a cylindrical surface. The glass beads are ballotini 
impact beads with a density of 2.45 gem -3 and diameters between 280 /mi to 300 ^im. 

One end of the tube is attached rigidly to a dc-motor, which rotates the tube around its 
cylinder axis. The motor speed is controlled with an accuracy of 10~ 4 . The control unit allows 
rotation frequencies / ro t from Hz to 50 Hz in steps of Hz with a maximal acceleration of 
1047 s -2 . The axial sand distribution in the tube is measured via light transmission, with a 
neon light placed behind the tube. A CCD camera, connected to a frame grabber, detects 
the intensity of the transmitted light. Each frame consists of a narrow strip (512 x 5 pixels) 
along the tube displaying high sand concentrations as dark stripes in a bright surrounding 
(Fig.lHb) and (c)). 

The experimental protocol consists of a stepwise increase or decrease of the rotation fre- 
quency from 2.4 Hz to 4.0 Hz followed by a subsequent waiting period of 5min, after which 
transients have disappeared. No patterns are visible for high frequencies above 3.9 Hz and 
frequencies below 2.9 Hz (Fig.EJa) and (b)). At intermediate frequencies patterns appear, 
depending on the history of the system. Thus we display the intensities for increasing and 
decreasing rotation frequencies in part (a) and (b) of Fig. [21 

Alternatively we characterize the pattern by a discrete Fourier decomposition of the lines 
(Fig. |2J. As soon as the rings appear, well defined peaks show up in the power spectrum. In 
the case of decreasing rotation frequencies the first peak is located at gi ow = 1.49 ±0.05 cm -1 , 
which corresponds to a ring spacing of 4.23 ± 0.15 cm. For frequencies below 3.5 Hz the first 
peak almost vanishes and a peak at (/high = 3.19 ± 0.05 cm -1 becomes dominant, due to a 
splitting of the rings. The pattern returns to <7i ow for frequencies between 3.4 Hz and 3.25 Hz. 
Below 3.25 Hz the particles start to accumulate at the right hand side of the tube. Whether 
the system chose the right or the left side could not be predicted in advance. For rotation 
frequencies below 2.9 Hz the sand distributed again homogeneously in the tube. Increasing 
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Fig. 2 - Different patterns at quasistatic changes of the (a) increasing and (b) decreasing rotation 
frequencies. Each line corresponds to one measurement of the sand distribution along the axes of 
rotation, which is an average of 20 images taken within 0.8 s. These intensities are normalized by 
the intensities obtained at zero rotation frequency, in order to prevent errors due to inhomogeneous 
illumination. The lines are plotted with an offset proportional to the rotation frequency to show the 
evolution of the system. The frequency was swept between 2.4 Hz and 4.0 Hz. (c) and (d) show the 
corresponding power spectra for increasing and decreasing rotation frequencies. Each line of (a) and 
(b) was transformed and plotted by the same offset as before. 

the rotation frequency leads to a different scenario as shown in Fig. EJa) and|2c). There the 
characteristic ring spacing is almost constant at q\ ow . No accumulation at a tube end was 
observed even after waiting periods of 30 min. 

To get a more precise picture of the transition points, we take as an order parameter 
the root mean square A rms of the measured intensities along the axes of rotation (Fig. |3J. 
We observe two hysteresis loops, which are related to the transitions between structured and 
unstructured states. These hysteresis loops indicate that the transitions can be considered as 
subcritical bifurcations. 

In order to characterize the bifurcation around 3.7 Hz further we now concentrate on the 
temporal evolution of the instability. We assume 

d t A = a A (1) 

where A describes the amplitude for a given wave number of the pattern, and a the growth 
rate. To measure a as a function of the driving frequency we start at a rotation frequency, 
where the granular medium forms a homogeneous layer (/ ro t = 4 Hz). Then we decelerate the 
motor suddenly to a frequency where the ring pattern is stable (/ ro t < /rot,c)- Due to this 
sudden change, the pattern will grow in time, as described by equation to lowest order. 
Figure 0] shows the growth of the amplitude in time after a jump from 4.0 Hz to 3.583 Hz 
(Fig-EIb)), and to 3.550 Hz (Fig. El^a)). The amplitude in this plot was obtained from the 
power spectrum by selecting the mode q\ ow = 1.49 cm -1 , which is the first stable one if the 
rotation frequency is decreased quasistatically (see Fig. Efd)). The solid line is the fit of an 
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Fig. 3 - Amplitude of the pattern as a function of the rotation frequency. 



exponential function to the data. At higher amplitudes the growth is no longer exponential. 
In fact, the amplitudes decays after 100s (Fig. Eta)). This non-monotonic effect is caused by 
a rearrangement of the pattern towards a higher wave number. 

This quenching of the system was performed for final rotation frequencies between 3.4 Hz 
and 3.9 Hz. Above 3.65 Hz no growth of the pattern was observed within a measurement 
time of 5 minutes. Below 3.4 Hz it was not possible to measure a growth rate for q\ ow . For 
frequencies between 3.633 Hz and 3.567 Hz the final stable mode was <7i w At lower frequencies 




Fig. 4 - Growth of the amplitude of the pattern as a function of time for two different rotation 
frequencies. The solid line represents the fit of an exponential function to the data. The insets show 
the regions of exponential growth on a semi-logarithmic scale. 
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this mode first grew in time and decayed afterwards, being replaced by the higher stable mode 
9high- Nevertheless, it was also possible to extract a growth rate for these frequencies, but 
with less accuracy, as indicated by the error bars in Fig. [5j 

All measured growth rates are plotted as a function of the rotation frequency in Fig. [5] For 
/ ro t between 3.55 Hz and 3.65 Hz the growth rate of the mode q\ ow increases roughly linearly 
with decreasing frequency, and saturates around 3.4 Hz, due to the different stable mode 9high, 
as seen in Fig. EJa). To model the increasing growth rate as well as the saturation, we fitted 
a parabola to the data. This parabola was then used to extrapolate the data for a — > 0. The 
corresponding critical rotation frequency is / ro t,c = 3.641 ± 0.035 Hz. 
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Fig. 5 - Growth rate as a function of driving frequency / ro t. 



Discussion and conclusion. - In order to estimate an upper boundary for the critical 
rotation frequency / ro t,c we analyze the motion of a single particle in the laboratory frame, 
neglecting the interaction with the neighboring particles. The calculation estimates the fre- 
quency for the detachment of a particle from the wall, which is a necessary, but not a suffi- 
cient criterion for pattern formation. A complete understanding of the instability cannot be 
achieved by this simple calculation. At high rotation frequencies the suspension is in rigid 
body rotation, and the particle sticks to the inner wall of the tube. Below a critical frequency 
the particle detaches from the wall and falls on a parabola-like trajectory described by the 
equation of motion 

mx = F G , res + F S + F cb ■ (2) 
-Fb,res represents the force of gravity, reduced by buoyancy: 

Fg, ics = ^nr 3 (pp - p f ) g . (3) 

r is the radius of the particle, p p and pf are the density of the particle and the fluid, respectively. 
The drag of the fluid on the particle is modelled by a modified Stokes formula 

Fs = —QTtnr (v p — Vf) - A , (4) 

where n denotes the dynamic viscosity of the fluid. v p and Vf are the velocity of the particle 
and the fluid, respectively, at the position of the particle. For the fluid velocity Vf we assume 
that the fluid is in rigid body rotation. Since the particle velocity is approximately equal 
to the fluid velocity, the Reynolds number is close to zero. A is the correction of the Stokes 
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Fig. 6 - (a) Trajectories of a particle falling in a viscous fluid for different rotation frequencies. The 
fluid is assumed to perform a rigid body rotation, (b) Maximal distance to the tube wall of a particle 
falling in a viscous fluid as a function of rotation frequency. The dashed line is drawn at the particle 
diameter. The inset is a close-up of the critical region. 



formula for particles close to the wall, derived by Brenner for vanishing Reynolds number |15j : 

2 sinh(2n + l)a + (2n + 1) sinh 2a 



A = — sinh a 7- ^7- ^ — — 

3 ^ (2n- 1 (2n + 3) 



4sinh 2 (n+ \)a - (2n + l) 2 sinh 2 a 



1 



(5) 



where the particle radius dependence which is taken care of by a is defined as 



a = ^h- 1 f'*±^) . (6) 



d is the distance of the sphere's surface to the wall. 

The third force is caused by the radial pressure distribution of the fluid in rigid body 
rotation: 

F ch = --TTr 3 pfUj 2 x , (7) 

the centrifugal buoyancy 

While the detachment point can be calculated analytically , the path of the particle is 
obtained by numerical integration of the equation of motion (J2J. Note that the particle, at 
the detachment point, has a distance oi R = R — r — 6 from the center of the tube. R is the 
inner radius and 6 the surface roughness of the tube. It is necessary to include the surface 
roughness, since the wall correction JSJ becomes singular for vanishing distance d. We have 
chosen 6—1 /im, which is a typical value for cast processes. 

Starting from the detachment point, we integrate the equation of motion J3J) until the 
particles touched the wall of the tube again. Three different trajectories are shown in Fig.EJa). 
From these trajectories we calculate the distance of the particle from the tube wall as a function 
of time. Each curve has a unique maximum d ma x, which depends on the rotation frequency of 
the tube. We determined numerically this maximal distance for different rotation frequencies 
(Fig-HHb)). For the rearrangement of the particles, which initially tend to form a close packed 
monolayer, into separated rings, some individual particles have to jump over their neighbors. 
In order to do this they have to achieve at least a distance d max equal to one particle diameter. 
If the distance is less than the diameter, particles will not be able to pass their neighbors. 
Therefore the critical frequency is found at the distance c? max equal to one particle diameter 
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2r. By this we obtain a critical rotation frequency of 3.642 ± 0.060 Hz, which is in agreement 
with the experimental value of 3.641 ± 0.035 Hz. 

A detailed theoretical investigation of the interaction between sand and water leading to 
the ring pattern remains to be done. A two-fluid approach might be difficult, since the primary 
flow is not trivial 18 . A molecular-dynamics simulation taking into account both the 
trajectories of the grains and the fluid flow, seems to be more promising. 

* * * 
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